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We present an update of the MILC studies of the physics of light pseudoscalars using improved 
staggered fermions. New runs at lighter quark mass, as well as increased statistics for older 
sets, are enabling us to improve the results for decay constants in full QCD. In addition, we 
have analyzed quenched runs at two different lattice spacings. This makes possible a test of the 
applicability of staggered chiral perturbation theory in a different context. 
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Since our previously published work [|I|] on the decay constants and masses of the 7i-K sys- 
tem, we have continued to generate and analyze lattices with Nf = 3 flavors of dynamical quarks, 
increasing the statistics and moving to lighter u,d, and (separately) lighter s masses. We use im- 
proved staggered fermions [Q]. With the same action, we have also investigated these quantities in 
the quenched approximation. Assuming the chiral fits are good and the continuum extrapolations 
are accurate, such calculations make possible: 

• A sensitive check of algorithms and methods — including the \/Det trick for dynamical 
staggered quarks — by comparing full QCD f n to the well-determined experimental value. 

• A precise extraction of the CKM matrix element V m from fx or fx/ fa, competitive with the 
world-average from alternative methods. 

• A determination of the light quark masses and their ratios with high lattice precision. 

• A test of the applicability of staggered chiral perturbation theory (SZPT) [pi ||] for describing 
lattice data — both Nf = 3 and quenched. 

• An extraction of the physical coefficients of the &(p 4 ) chiral Lagrangian. 

• A determination of the extra, unphysical parameters that enter SZPT. This is important be- 
cause these parameters also appear in staggered chiral forms for other physical quantities, 
e.g., heavy-light decay constants and form factors Fixing the parameters from the light- 
light system reduces the systematic errors in heavy-light computations [Q]. 

In our computations, we analyze two lattice spacings: a«0.12 fm ("coarse") and a «0.086 fm 
("fine"). We call the valence quark masses are m x and m y ; the dynamical quark masses (for Nf = 3) 
are m' u = m' d = m', and m' s . Here primes indicate simulation values; the corresponding masses 
without primes are the physical values. Table [T] shows the lattice parameters used. 

The relative lattice scale is determined using the length r\ [Q] from the static quark potential. 
We reduce statistical fluctuations in r\/a by fitting to a smooth function. The absolute lattice scale 
is obtained from the T 2S- IS splitting. Following the continuum extrapolation in Ref. []8|], but using 
updated HPQCD results [fl, rather than those in [fill], we obtain r\ = 0.318(7) fm. 



am q 1 am s 


10/g 2 


size 
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number 


a (fm) 


0.03 / 0.05 


6.81 


20 3 x 64 


(2.4 fm) 3 


564 


0.120 


0.02 / 0.05 


6.79 


20 3 x 64 


(2.4 fm) 3 


484 


0.120 


0.01/0.05 


6.76 


20 3 x 64 


(2.4 fm) 3 
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0.121 


0.007 / 0.05 


6.76 


20 3 x 64 


(2.4 fm) 3 


493 


0.121 


0.005 / 0.05 


6.76 


24 3 x 64 


(2.9 fm) 3 


400 (197) 


0.120 


0.03 / 0.03 


6.79 
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(2.4 fm) 3 


350 


0.120 


0.01 / 0.03 


6.75 


20 3 x 64 


(2.4 fm) 3 


349 


0.120 


quenched 


8.00 


20 3 x 64 


(2.4 fm) 3 


408 


0.119 


0.0124/0.031 


7.11 


28 3 x 96 


(2.4 fm) 3 


531 


0.0863 


0.0062/0.031 


7.09 


28 3 x 96 


(2.4 fm) 3 


583 


0.0861 


0.0031 / 0.031 


7.08 


40 3 x 96 


(3.4 fm) 3 


210 


0.0859 


quenched 


8.40 


28 3 x 96 


(2.4 fm) 3 


396 


0.0853 



Table 1: Lattice parameters. Runs and numbers of configuration in normal font were included in Ref. [|ljj 
those in bold font are new. The lattice sets above the double line are "coarse;" those below are "fine." 
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As in Ref. [{[]] we fit the partially quenched lattice data to SZPT forms [Eh. The addition of 
new runs allows us (sometimes forces us) to change some details of the fits. To determine LO and 
NLO chiral parameters we fit only to the low quark-mass region. The cut on valence quark mass 
is the same as before: am x ± am y < 0.021 ~ 0.5am s (coarse) and am x + am y < 0.017 « 0.6am s 
(fine). We now have enough data to cut on sea-quark masses, too: We remove the 0.03/0.05, 
0.02/0.05, and 0.03/0.03 sets for this fit. Because statistical errors are so small, (0.1% to 0.4% for 
decay constants, 0.1% to 0.7% for squared meson masses), we still require NNLO analytic terms 
in addition to complete NLO forms to get good fits. Such joint fits to decay constants and masses, 
including both coarse and fine lattices, have 26 free parameters: 

• 2 LO parameters: / and pL (decay constant and condensate at tree-level). 

• 8 NLO parameters: 4 physical and 2 taste-violating analytic terms, 2 taste-violating hairpins. 

• 10 physical, NNLO analytic parameters. 

• 6 tightly constrained parameters (prior width = 0.025): give variation of 2 LO and 4 NLO 
physical parameters with lattice spacing. 

For interpolation around m s , we must include higher quark masses. Once the LO and NLO 
parameters are determined, we fix them (up to statistical errors) and fit to all sea mass sets and wider 
ranges of valence masses. For central values we choose the range am x + am y < 0.055 l.4m s 
coarse, and am x + am y < 0.0353 rj 1.3m s fine. This fit is called "Fit C." For systematic error tests, 
the range is widened to am x + am y < 0. 10 2.5m s coarse, and am x + am y < 0.062 « 2.2m s fine. 
With either of these choices, we need to add in the NNNLO analytic terms (18 parameters, cubic 
in quark masses for f K and M\j (m x + m y )) to get good fits. 

With our old data set, m' s only changed with a, and m' s was usually significantly larger than m', 
so the sea quark mass dependence and a dependence were difficult to disentangle. The new data, 
which includes coarse lattices with am' s = 0.03 in addition to the previous value am' s = 0.05, gives 
better control of the sea quark mass dependence and smaller a dependence of the LO, NLO, and 
NNLO parameters. Including the NNNLO terms and the a dependence of the NNLO terms gives 
28 parameters more than the low-mass fits described above, for a total of 56 parameters. Twelve of 
these (LO, LO a dependence, and NLO parameters) are tightly constrained from the low-mass fits. 

Figure [j] shows Fit C results for f n and fg. This fit has # 2 /d.o.f. =0.99 for 556 degrees 
of freedom (confidence level CL=0.59). In each plot, the maroon line is the "prediction" for the 
0.0031/. 031 fine run based on a second fit that leaves out that run; it should be compared with the 
solid black line that comes from Fit C. Since the time of the conference, we have accumulated about 
25% more 0.003 1/.031 configurations, and the effect of removing or including this run in the fit has 
decreased. Given that the CL decreases when the run is removed, we no longer see any reason to 
consider omitting the run. We note that the 0.003 1/.031 run is still only about half finished, so there 
will probably be further noticeable shifts. In the plot, two "experimental" points (shifted slightly 
to the left for clarity) are shown. Both points are based on the measured leptonic (K — > Iv) rate, but 
* assumes the PDG value V us = 0.2200(26) [jfl]]; while x assumes the results of recent experiments 
V us = 0.2262(23) [|l2|]. Both these values of V us come from experimental determinations of the 
semileptonic (K — » niv) rate and non-lattice theory for form factors. 

Our preliminary results for decay constants are: 

f n = 128.1 ±0.5±2.8MeV, f K = 153.5 ±0.5 ±2.9 MeV, f K /f K = 1.198(3)(± 1 |) , 
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Figure 1: Left: Comparison of Fit C to partially quenched f K data. Extrapolation to continuum, setting 
m' s = m s , and going to full QCD gives the red line. Red + shows the final result after extrapolation m x ,m y ^> 
m. The maroon line (just barely visible above the black line) is the prediction for the 0.0031/0.031 fine run 
(black squares) based on the other data. Right: Same as left, but for fx- The short green continuation of the 
red line keeps light sea masses fixed at the average physical value m and extrapolates m x —>m u , 

where the errors are from statistics and lattice systematics. These results are consistent with our 
previous answers [fi]], with slightly smaller errors. The current Nf = 3 results for quark masses are 
little changed from those in Refs. [jj], |13| ]. 

We extract V us from our ft:/ fit result. This is probably safer than using fx itself, because the 
ratio is largely free of scale errors. We obtain \V US \ = 0.2242(^3}), which is consistent with world- 
average values, with comparably sized errors. From Fig. [T], one can deduce that using fg alone 
would result in a somewhat higher value of V us . The difference comes from the fact that our f n 
result is slightly low compared with experiment, although consistent within errors. Runs planned 
for the near future, as well as those now in progress, should allow a further reduction in the errors. 

We now turn to the quenched data. We fit to the same valence mass range as the Nf = 3 Fit C. 
Again, terms through NNNLO are needed; joint fits to decay constants and masses, including both 
coarse and fine lattices, have 34 free parameters: 



• 3 LO parameters: /, jj,, and the quenched chiral parameter 8 Q14|]. We consider 8 to be "LO" 
because its effects are not suppressed by powers of quark mass. 

• 7 NLO parameters: 2 physical and 2 taste-violating analytic terms, 2 taste-violating hairpins, 
and the quenched chiral parameter a [p4|]. 

• 4 physical, NNLO analytic parameters. 

• 4 physical, NNNLO analytic parameters. 

• 14 tightly constrained parameters (prior width = 0.04): give variation of 2 LO, 4 NLO, 4 
NNLO, and 4 NNNLO physical parameters with lattice spacing. 

• 2 parameters to allow the r\ scale on the coarse and fine lattices to vary within la. 

The parameter 8 multiplies a function of the taste-singlet mass, which is large ( > 500 MeV) on 
coarse lattices because of taste splitting. The coarse lattices are therefore insensitive to quenched 
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Figure 2: Left: Comparison of Fit Q to quenched / w data. Continuum extrapolation gives the red curve. 
(ra„ + mrf) phys is from the Nf = 3 analysis. Right: Same as left, but for quenched fx. 

chiral logs, and 8 is poorly determined from coarse-lattice fits: 0.02 <<5< 0.18, depending on 
details of the fit and the higher order terms included. Currently, our preferred approach is to obtain 
8 from fits to the fine lattices alone, and hold it (as well as the fine lattice scale) fixed in the joint 
fits. The resulting joint fit, called Fit Q, has 32 parameters and # 2 /d.o.f. = 0.96 for 98 degrees of 
freedom; CL=0.59. Comparison of this fit with the quenched f n and f K data is shown in Fig. |2| 

The continuum-extrapolated curve (red line) blows up at low quark mass for f K but not for f n : 
the latter does not have a 1-loop quenched chiral log in the continuum. Taking into account SZPT, 
f n increases significantly with the a — > extrapolation. This makes makes our quenched f % differ 
from experiment (or our Nf = 3 results) by an even larger amount than was reported in Ref. [10] 
from coarse data alone. In fact, we find /^ uench //^ /_3 ~ 1.28 (with the r\ scale). Note, however, 
that the difference between the raw coarse and fine data for f n is small, so it is not yet clear how 
seriously we should take the continuum extrapolation — data at smaller a is needed. This is an even 
more important issue for quenched fy, where the continuum result blows up in the chiral limit, but 
the lattice data is smoothed out by SZPT effects. (This is an example of the lack of commutativity 



of chiral and continuum limits in SZPT for infrared sensitive quantities [15].) 

Values of the taste-violating hairpin parameters in the quenched analysis are comparable to 
those for Nf = 3. We find: 



-0.29(1)(4) 
-0.12(2)(+i 1, 



r\a 2 {8' A )^ ch = -0.20(2) (+^) r\a 2 {8' A ) N f=^ ■■ 

rja 2 (^)q uench = 0.09(5)(+^) r 2 a 2 (8{,) N f= 3 - s 

In both quenched and = 3, 8^ is poorly determined and consistent with 0. For the quenched 
parameter 8 we are finding: 8 = 0.09(1)(5). This result is consistent with most quenched evalua- 
tions [16], which get 8 « 0.1, but not that of the Kentucky group [17], who obtain 8 = 0.24(3)(4). 
(Results in Ref. [ |l8| ] also tend toward higher values.) Assuming the analysis in Ref. [|l7|], this 
difference makes sense. Despite the low quark masses in the current calculation, the large taste- 
violations in the taste-singlet sector puts us in a larger region of meson mass, where the effective 8 
from chiral fits is smaller. Modulo this issue, which needs further study, it seems that SZPT works 
reasonably well in the quenched theory, as it does in the Nf = 3 case. 
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